Theoretical study of electron states in Au chains on NiAl(llO) 
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We have carried out a density functional study of unoccupied, resonance states in a single Au atom, 
dimers, a trimer and infinite Au chains on the NiAl(llO) surface. Two inequivalent orientations of 
the ad-chains with substantially different interatomic distances were considered. From the study 
of the evolution of the electron states in an Au chain from being isolated to adsorbed, we find 
that the resonance states derive from the 6s states of the Au atoms, which hybridize strongly with 
the substrate states and develop a p-like polarization. The calculated resonance states and LDOS 
images were analyzed in a simple tight-binding, resonance model. This model clarifies (1) the 
physics of direct and substrate-mediated adatom-adatom interactions and (2) the physics behind 
the enhancements of the LDOS at the ends of the adatom chains, and (3) the physical meaning 
of the " particle- in-box" model used in the analysis of observed resonance states. The calculated 
effective mass and band bottom energy are in good agreement with experimental data obtained 
from scanning tunnelling spectroscopy. 

PACS numbers: 73.20.Hb, 73.21. Hb, 68.37.Ef 
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I. INTRODUCTION 

A most interesting development in nanoscience has 
been the possibility to tailor and analyze the electronic 
properties of single nanostructures of adatoms on metal 
surfaces, This development has been made possible by 
the unique capabilities of the scanning tunnelling mi- 
croscope to image and manipulate single adatoms and 
to probe their electronic structure by tunnelling spec- 
troscopy. A pioneering example is provided by the quan- 
tum corral experiments, in which atomic manipulations 
were used to assemble closed atomic structures, which in 
turn act as two dimensional resonance cavities for sur- 
face state electrons on a metal surface^. Another most 
interesting example is provided by the unoccupied, elec- 
tron resonance states in Au adatom and addimers^ and 
longer Au ad-chains^*^ assembled on a NiAl(llO) sur- 
face, which were revealed by imaging the differential con- 
ductance for various biases. This class of resonance states 
was also identified in linear Cu chains on Cu(lll)& so 
their existence is not limited to the NiAl(llO) surface. 

The appearance of such resonance states in Au adatom 
structures was rather surprising and raises several ques- 
tions about their physical nature. In the original experi- 
ments*^, the differential conductance spectra and images 
of linear Au ad-chains were analyzed using the quan- 
tum states of a simple, one-dimensional "particle-in-box" 
(PIB) model. The observed dispersion of the resonance 
states in the Au chain was found to be free-electron like 
with an effective mass being about half of the bare mass 
and the observed dependence of the low energy peak posi- 
tion on the ad-chain length followed closely the prediction 
of the PIB model. However, the origin and the character 
of these resonance states were not clarified in this model 
and the physical meaning of the one-dimensional "box" 
is unclear in view of the discrete atomic structure of the 
ad-chains. Thus there is a need to gain a more detailed 
understanding of these resonance states from calculations 



and modelling of their electronic structure. 

The unoccupied, resonance states of the Au adatom 
structures are excited states that are probed by adding 
an extra electron. In principle, they are not directly de- 
scribed by a ground state theory such as density func- 
tional theory. However, unoccupied, Kohn-Sham states 
provide a useful, zero-order approximation of such ex- 
cited states^ as demonstrated by density functional calcu- 
lations of a single Au adatom and addimers with various 
interatomic separations. These calculations showed that 
the observed resonance structures in the STS were repro- 
duced in the calculated local density of states (LDOS) at 
the tip apes^. These resonances were formed in an en- 
ergy region with a depletion of LDOS and had a mixed 
s and p character. 

A first attempt to understand the nature of the reso- 
nance states of the Au ad-chains was the recent density 
functional study of isolated, finite Au chains by Mills and 
coworkers 10 «. They argued that the observed resonance 
states derived from unoccupied states with ir character 
of the isolated chains. This argument resolved the issue 
that the observed lowest energy state with no nodes along 
the chain is unoccupied in contrast to the states of the iso- 
lated chains with a character. Furthermore, they intro- 
duced a three-dimensional "particle-in-a-cylinder" model 
to rationalize all STM observations. However, this study 
did not address the effects of the interaction of the states 
in the Au chains with the substrate states. 

In this paper, we present a density functional study 
of the resonance states of some linear Au adatom struc- 
tures adsorbed along two inequivalent directions on the 
NiAl(llO) surface. The focus is primarily on the infi- 
nite Au ad-chains of which one has the same orienta- 
tion along the [001] direction as an assembled chaini£, 
whereas the other infinite, Au ad-chain is oriented along 
the [110] direction. The substantially larger interatomic 
distance of the latter ad-chain than the former ad-chain 
illustrates the effects of interatomic interactions on the 
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electronic states. We present also results for a single 
adatom, two different ad-dimers and an ad-trimcr. The 
nature of the resonance states in the infinite ad-chains is 
revealed by studying their evolution from the states of the 
infinite, isolated Au chains with decreasing chain-surface 
distance. We also discuss the calculated resonance states 
and the tunnelling through the resonances states in terms 
of a tight-binding, resonance model. This model clarifies 
the physical meaning of the "particle-in-box" model in 
the analysis of the resonance states. 

The paper is organized in the following manner. In 
Section |nj we present some details of the density func- 
tional calculations of the adatom structures and their 
geometric and electronic structure are presented in Sec- 
tion Ifffl The nature of the resonance states in the Au 
ad-chains and their relation to the states of the isolated 
Au chains are discussed in Section lfV Al A simple, tight- 
binding resonance model for the states in the adatom 
structures is described in the Appendix and discussed in 
relation to calculated resonance energies and images in 
Section If VBI Our results are compared with experimen- 
tal data in Section where we also clarify the physical 
meaning of the "particle-in-box" model in terms of the 
tight-binding, resonance model. Finally, we give some 
concluding remarks in Section IVII 



II. COMPUTATIONAL METHOD 

We have considered several linear Au adatom struc- 
tures on the NiAl(llO) surface including single adatoms, 
two different addimers and infinite, adchains along along 
either the [OOf] direction or the [110] direction and an 
adtrimer along the [001] direction. The [001] and the 
[110] directions correspond to the T — Y and the T — X 
directions in the surface Brilliouin zone (SBZ), respec- 
tiveljiii so that the corresponding chains aligned along 
these two directions are henceforth referred to as the Y 
and X chain, respectively. The density functional calcu- 
lations of the single adatoms and the two addimers have 
already been briefly described in Ref. 0- 

The electronic and the geometric structure of the Au 
chains and the trimer on NiAl(llO) were studied us- 
ing density functional calculations that were carried out 
using the projector augmented wave method as imple- 
mented in the VASP codei2iiii4. The exchange and cor- 
relation effects were represented by the generalized gra- 
dient approximation— The two infinite, (periodic) chain 
structures were represented by an Au adatom on a NiAl 
slab in two different super cell geometries. To minimize 
the interactions between the chain and its periodic im- 
ages, we used an 1x4 and 6x1 surface unit cells for the Y 
and X chains, respectively. The results for the electronic 
states of the Y chain were found to be well-converged for 
9 substrate layers and this number of layers was also used 
for the X chain. The vacuum region of the super cell was 
6 layers. For the adtrimer we used a super cell with a 5x2 
surface unit cell and 8 layers of NiAl. The kinetic energy 



cut-off for the plane wave basis set was about 200 eV and 
the surface Brillouin zone were sampled by 80, 28, and 
24 k points for the Y, X ad-chains, and the adtrimer, 
respectively. The Au atom and the two outermost layers 
of Ni and Al atoms were fully relaxed until the residual 
forces were less than 0.05eV/A. 

To make contact with scanning tunnelling spectra, we 
have calculated the local density of states (LDOS) out- 
side the surface. This approach is based on the Tersoff- 
Hamann (TH) approximation for tunnelling in an STM 
junction 16 . In this approximation, the differential con- 
ductance, 4{r(V), as defined from the tunnelling current 
J as a function of sample bias V, is given by, 

^(V) ex p(r ,e = eV) (1) 
= ^|^(r )| 2 ^-^) (2) 

Here the local density of states (LDOS), p(r ,e — eV), 
at the position rg of the tip apex is expanded in wave 
functions Vv^o) °f Kohn-Sham states /i with energy 
of the sample in the absence of the tip. The continuum 
of one-electron states is mimicked in the slab calculations 
by a Gaussian broadening of the delta function in Eq. |5J 
In principle, the range of applicability for TH approxi- 
mation to the differential conductance is limited to small 
V. At larger V, one has in principle to take into account 
both the change of the electronic states by the electric 
field from the tip and the V dependence of the tunnelling 
barrier. As suggested by Langii, we have accounted for 
the latter effect by extending the wave functions into the 
vacuum region using a decay parameter set by the In- 
dependent vacuum barrier half-way from the surface to 
the tip apex. 



III. RESULTS 

In this section we will present our results for the LDOS 
of the various adatom structures in the progression of 
adatoms, ad-dimers, linear ad-trimers to infinite chains 
of Au atoms on the NiAl(llO) surface. The effects of the 
adatom-adatom distance on the resonance states are re- 
vealed by considering the alignment of the addimer and 
the ad-chain along the two inequivalent X and Y direc- 
tions of the NiAl(llO) surface as defined in Section ITT1 

We begin by presenting the results for the geometries 
of the various linear adatom structures studied in this 
work. Following the suggestion by the STM experiments, 
we started our structural optimization by placing the Au 
adatoms in the short Ni bridge sites. In the optimiza- 
tion for the ad-dimers and the ad-trimer, we find that 
the adatom-adatom interactions are weak compared to 
the adatom-substrate interactions so that all adatoms 
are located laterally within 0.03 A from the short Ni 
bridge positions. Thus the adatom-adatom distance is 
about 2.9 A and 4.1 A for the linear, adatom structures 
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FIG. 1: Local density of states (LDOS) of (a) a single Au 
adatom and (b) Au ad-dimers on NiAl(llO), and (c) Au ad- 
trimer. (b) Solid and dashed lines is the LDOS for the ad- 
dimer along the Y and X directions, respectively. The LDOS 
corresponds to a tip-surface distance of about 7 A and in a 
position either (a) on top of the adatom, (b) laterally dis- 
placed 1 A away from an adatom along the dimer axis, or (c) 
on top of an ed ge a datom. The measured peak energies as 
taken from Refs. l3l7l are also indicated by vertical bars 



along the Y and the X directions, respectively. How- 
ever, the structures with the short adatom-adatom dis- 
tance tend to increase slightly the adatom-surface dis- 
tance and relax slightly some surface Ni atom positions. 
The adatom-surface distance increases from 1.95 A for 
the single adatom to 2.05 A for the ad-chain along the 
Y direction. The bare NiAl(llO) surface was found to 
be rumpled with the Al rows being displaced 0.18 A 
farther out from the surface than the Ni rows in good 
agreement with earlier experimental and theoretical stud- 
iesim. The only substantial adatom-induced substrate 
relaxation was found for the linear structures along the Y 
direction and involved a downward relaxation by about 
0.1 A of the Ni atoms that are coordinated to two Au 
adatoms. 

The results for the LDOS of the bare NiAl(llO) sur- 
face, the single adatom and the two ad-dimers are shown 
in Figs.Q] As already reported and discussed in Ref . 0, 
the LDOS of the single Au adatom exhibits a single res- 
onance peak in the unoccupied LDOS at about 1.71 eV, 
which is split into a resonance doublet in the LDOS for 
the two ad-dimers. The result for the LDOS of the clean 
surface showed a depletion in the structure-less LDOS in 
this energy range and exhibited a sharp onset around 2.5 
eV above the Fermi energy. The physical origin of these 
resonance states will be discussed further in Section flV Al 
From the character of the LDOS images at the peak en- 
ergies, the low-lying state in the resonance doublet of the 
ad-dimer was shown to be a symmetric combination and 
an anti-symmetric combination of the resonance states of 
the single adatoms in accordance with a simple, two-state 



FIG. 2: Local density of states images of resonance states for 
the Au trimer on NiAl(llO). (a) First state at 1.0 eV, (b) 
second state at 1.6 eV, and (c) third state at 2.0 eV . The 
origin is located at the center atom of the ad-trimer. Same 
tip-surface distance as in Fig. 



model. The larger resonance splitting for the ad-dimer 
along the Y direction than along the X direction is sim- 
ply caused by the stronger adatom-adatom interaction 
for the shorter ad-dimer than for the longer ad-dimer. 

In the case of the linear ad-trimer, the interactions 
among the three adatom resonances give rise to a reso- 
nance triplet, as shown in Fig.^c). The maximum reso- 
nance energy splitting is about 1.0 eV and is larger than 
the energy splitting of 0.9 eV for the Y ad-dimer. The 
characters of the resonance states are also revealed by 
LDOS images at the resonance peak energies, as shown 
in Fig. |21 Note that the image for the lowest lying reso- 
nance state has no nodal planes, whereas the other states 
has an increasing number of nodal planes with increas- 
ing resonance energy. The resonance energy positions 
and characters of these states and the strong intensity 
of the end lobes in the LDOS image of the highest lying 
resonance state will be discussed in Section ILV Bl 

In the case of longer ad-chains, the interactions among 
the increasing number of adatom resonance states should 
eventually give rise to a band of resonance states for the 
infinite ad-chains. The formation of such a band is shown 
in the LDOS in Fig.[31for the infinite ad-chains along the 
Y and X directions. The stronger adatom-adatom inter- 
actions for the Y ad-chain than for the X ad-chain results 
in a larger bandwidth width for the Y ad-chain than for 
the X ad-chain. The dispersions of the bands of reso- 
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FIG. 3: Local density of states of the Au chains on NiAl(llO) 
along the (a) Y and (b) X directions. The tip is above an Au 
atom at the same distance as in Fig. 




Reduced wave vector 



FIG. 4: Wave-vector-resolved local density of states of the 
Au chains on NiAl(llO) along the (a) Y and (b) X direc- 
tions. The dashed lines are dispersions as obtained from a 
nearest neighboring tight binding model, whereas the solid 
lines are the corresponding free particle dispersions with the 
same effective mass. 

nance states are revealed by the resolution of the LDOS 
over wave vectors in the one-dimensional Brilluion zones 
(LBZ) of the two ad-chains. In Fig. we show contour 
plots of these wave- vector resolved LDOS. These results 
show that the LDOS exhibit a resonance structure that 
disperses with the wave vector up to the middle of LBZ 
where the dispersion levels off and the resonance becomes 
ill-defined. For small wave vectors the dispersion is free- 
electron like with an effective mass of 0.65 and 0.91 m e 
for the Y and X ad-chains, respectively. 



IV. DISCUSSION 

In this section we will discuss the physical origin of 
the resonance states in the Au adatom structures on 
NiAl(llO). This discussion concerns primarily the for- 



mation of the resonance states in the one-dimensional 
Au ad-chains. The energy positions and the LDOS im- 
ages of the resonance states of the adatom structures will 
be discussed in a simple, tight-binding, resonance model. 
This model enables us also to discuss the LDOS images 
of ad-chains with arbitrary lengths. 



A. Nature of the resonance states in the Au 
ad-chains 

The existence of resonance states in Au adatom struc- 
tures on the NiAl(llO) surfaces raises several questions 
concerning the relation of these states to the states of 
the isolated adatom structures and the role played by 
the electronic states of the substrate. These questions 
are addressed here by scrutinizing the resonance states of 
the infinite ad-chains. We begin by discussing the states 
of the isolated Au chains. 

In Figs. Ufa) and (b), we show the calculated band 
structures for the isolated Au chains with same geometry 
as the ad-chains along the Y and the X directions. These 
results are easily understood in terms of the electronic 
structure of an isolated Au atom as obtained from spin- 
unpolarized density functional calculations. The five- 
fold, spatially degenerate 5d states of an Au atom are 
fully occupied with an energy 1.8 eV below the energy 
of the singly occupied 6s state. The three-fold, spatially 
degenerate 6p states are unoccupied and located about 5 
eV above the 6s state. When forming an infinite, chain 
these atomic states overlap and form bands. In the case 
of the X chain, the interatomic distance is about 4.1 A 
and the atomic states form narrow and non-overlapping 
bands with predominantly atomic character as shown in 
Fig. |3{b). Only the 6p states show an appreciable band 
width because they are close in energy to the vacuum 
level so that the corresponding wave functions are ex- 
tended and have a large overlap. For the Y chain with a 
shorter interatomic distance of about 2.9 A, the disper- 
sion of the 5d and 6s states increase substantially com- 
pared to the X chain and develop a mixed character. 
The behavior of these band structures does not provide 
a simple answer to the question which atomic states are 
involved in the formation of the unoccupied resonance 
states in the ad-chain. The 6p bands are unoccupied but 
are far away from the Fermi level and the 6s band is 
half-occupied. 

To understand how the resonance states of the Au ad- 
chains are related to the states of the isolated Au chains, 
we have calculated the evolution of the LDOS and the 
partial DOS for an Au atom in the Y ad-chain at the 
r point of the LBZ as a function of the outward rigid 
displacement Ad cs of the chain from its equilibrium po- 
sition. In Fig. El we show the calculated LDOS at a fixed 
distance of 3 A outside the chain as a function of Ad cs . 
At the largest Ad cs of 2.1 A, the chain-substrate cou- 
pling is weak and the LDOS have a prominent peak that 
is close in energy to the occupied states of s character of 
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FIG. 5: Band structure of the isolated Au chains with same 
interatomic distances as the ad-chains along the (a) Y and 
(b) X directions. Bands with predominantly s, p, and d char- 
acter are marked by solid circles, diamonds and crosses, re- 
spectively. 
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FIG. 6: Evolution of local density of states of the infinite, Au 
chain on NiAl(llO) along the Y direction with outward dis- 
placement Adas of the ad-chain from its equilibrium distance. 
The tip-ad-chain distance is kept at a fixed distance of 3 A 
above an Au atom. 

the isolated chain. For decreasing d cs , this state of domi- 
nant s character of the chain shifts upward in energy and 
broadens into the unoccupied resonance state at the equi- 
librium position of the chain. Thus there is no indication 
in these results for the LDOS that this resonance state 
evolve from a 6p z state of the Au chain that is broaden 
into a resonance and shifted down continuously in energy 
with decreasing Ad cs . However, the result that the oc- 
cupied 6s state of the chain at the Y point turns into a 
narrow unoccupied resonance upon adsorption appears to 
be indicate that the ad-chain is not charge neutral. This 
conflicting result is resolved by scrutinizing the evolution 
of the partial DOS. 

In Fig. [3 we show the evolution of the partial DOS 
for an Au atom with s, p z , and d z 2 characters at the 




-5 -4 -3 -2 -1 1 2 3 4 5 

e-£ F (eV) 

FIG. 7: Partial wave decomposition of the density of states 
for an Au atom in the chain on NiAl(llO) along the Y direc- 
tion for different outward displacements, Ad cs , of the ad-chain 
from its equilibrium ad-chain-surface distance, as indicated in 
the figure. 



r point in the chain BZ as a function of Ad cs . At the 
largest Ad cs of 2.1A, the chain is essentially decoupled 
from the substrate and the partial DOS exhibit peaks at 
energies close to the energies of the states of the isolated 
Au chain. The state at about 0.5 eV below the Fermi 
level has s character and broadens appreciably with de- 
creasing Ad cs and develops both a resonance structure 
in the unoccupied states and resonance structures in the 
occupied states well below the Fermi level. Thus the for- 
mation of both occupied and unoccupied states shows 
that the formation of the resonance with s character in 
the unoccupied states is not in conflict of the ad-chain 
being neutral. The partial DOS shows also that the res- 
onance state has a strong p z character at the equilibrium 
position but no significant d z 2 character. The develop- 
ment of the p z character of the resonance state with de- 
creasing Ad cs is not correlated with the 6p z state of the 
isolated, chain but involves rather a polarization. This 
latter state broadens substantially into a wide band upon 
adsorption with decreasing Ad cs and does not shift con- 
tinuously down in energy. 

The calculated band structure of the bare NiAl(llO) 
surface shows that there are projected band gaps in a 
energy region where the resonance state of the ad-chain X 
at the r point is formed. In Fig.[5| we show the calculated 
electronic states of a slab with 18 substrate layers along 
the Y-X and Y-Y directions in the surface BZ (SBZ) of 
NiAl(llO). Only states that are even with respect to the 
symmetry planes spanned by these directions in the SBZ 
and the surface normal are shown in Fig. El The T point 
in the LBZ of the ad-chain X corresponds to the line of k 
points along the Y — X direction in the SBZ. Along this 
line there is a projected band gap in an energy region 
that embraces the resonance energy about half way to 
the LBZ boundary which show that this resonance state 
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FIG. 8: Projected band structure of the isolated NiAl(llO) 
surface as a function of the reduced wave vector along the 
Y and X directions. Only states with even character with 
respect to these directions has been displayed. The bands 
were calculated for an 18 layer slab. 
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FIG. 9: Calculated bulk density of states of NiAl(llO) (solid 
line) and Cu (dashed line). 



is indeed a resonance since it overlap with the bulk states. 

A depletion of density of states in the energy region 
of the Au adatom-induced resonance states is revealed in 
the bulk DOS. In Fig. EI we compare the calculated DOS 
for bulk NiAl and Cu. In the energy range of about -5 to 
-2 eV the DOS is dominated by the d states. In the case 
of NiAl, the transition metal atoms are surrounded by a 
simple metal atom resulting in a much narrower d band 
than for Cu. At larger energies the DOS are dominated 
by states with sp character. In the case of NiAl there is 
a depletion in the DOS in energy range from about to 
2.5 eV, whereas Cu shows no such depletion of the DOS. 
Note that this depletion is not a necessary condition for 
the formation of resonance state because similar states 
have been shown to exist in Cu adatom structures on 
Cu(lll)&. In the latter case the resonance state is formed 
in the prominent sp band gap of the Cu(lll) surface. 



B. Simple model for resonance states and images 

To gain a better understanding of the interactions 
among the resonances states and the LDOS images for 
finite chains, we have analyzed our results in a simple 
tight binding, resonance (TBR) model with an s-wave 
approximation for the tails of the resonance wave func- 
tions in the vacuum region. This model is defined in the 
Appendix. The TBR model explains the success of and 
justifies the "particle-in-the-box" analysis of the observed 
scanning tunnelling spectra for the Au chains. 

We begin by investigating to what extent the TBR 
model can describe the calculated resonance energies of 
the various adatom structures. One obvious choice of 
tight-binding parameters is based on the calculated en- 
ergies of the resonance doublet of the ad-dimers. In the 
case of the ad-dimer Y, these energies are 1.2 and 2.1 
eV, which are reproduced by an on-site energy eo = 1.65 
eV and a nearest-neighboring off-site energy t\ — -0.45 
eV. Note that in this model eo corresponds to the res- 
onance energy for a single adatom and eo is very close 
to the calculated value of 1.7 eV for the adatom. The 
predicted values for the resonance energies for the reso- 
nance triplet of the ad-trimer are then eo -Mi/v2 = 1.01, 
eo=1.65, and eo — ti/\/2 = 2.29 eV. These values are in 
good agreement with the calculated energies for the two 
lowest-lying states but for the highest-lying state the cal- 
culated energy is too high with about 0.3 eV. In the case 
of the ad-dimer X, the resonance energies suggests a on- 
site resonance energy, eo = 1.70 eV that is close to its 
value for the ad-dimer Y, whereas the off-site energy t± 
= -0.30 eV is smaller than its value for the ad-dimer Y. 

In our earlier study of the resonance doublets in ad- 
dimers with various interatomic distances, we argued 
that there are both direct and substrate-mediated con- 
tributions to the adatom-adatom resonance interaction 
described by t\ . The origin of these contributions is clar- 
ified by the result of the TBR model for t\ in Eq. I A. 61 
t\ = tf+tf has a direct contribution and a substrate- 
mediated contribution tf. if = Zt^ is given by the 
direct interaction between the 6s atomic orbitals being 
renormalized by their resonance strength Z. decays 
exponentially with the interatomic distance d whereas tf 
is long ranged and was found to dominate t\ for larger d. 
tf can also have a negative imaginary part which would 
enhance and diminish the broadening for the symmetric 
and anti-symmetric resonance states, respectively. How- 
ever, this effect is not discernable in the calculated LDOS 
for the ad-dimers. 

In the case of the infinite ad-chains, the sets of TB 
parameters derived from the ad-dimers give a poor de- 
scription of the dispersions of the resonance states for the 
Y and X ad-chains. The predicted energies e + 2ti = 0.75 
and 1.10 eV of the lowest energy state of the TB model 
for the Y and X ad-chains, respectively, are in agree- 
ment with the calculated energies of 0.69 and 1.23 eV for 
the corresponding ad-chains. However, these parameters 
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gives the wrong trend for the effective masses m* for the 
ad-chains, as defined by the near-parabolic behavior of 
the dispersion for small wave vectors. The TB model 
with these parameters gives m* — h 2 /(2a 2 ti) — 1.0 and 
0.75 m e for the ad-chain Y and X, respectively, whereas 
from the calculated dispersion one obtains the reverse 
trend: m* = 0.65 and 0.91 m e for the ad-chain X and 
Y, respectively. The breakdown of this model in de- 
scribing both ad-dimers and ad-chains is a consequence 
of substrate-mediated long-range interactions among the 
resonance states, as was demonstrated by the combined 
theoretical and experimental study in Ref. for the res- 
onance states of ad-dimcrs with various interatomic dis- 
tances. However, as shown in Fig. the near-parabolic 
dispersion of the resonance states for wave vectors up to 
half of the LBZ is described well by the effective TB pa- 
rameters e = 2.09 (1.83) and t Y = -0.7 (-0.25) eV for 
the Y (X) ad-chain. Note that for larger wave vectors, 
the strengths of the resonance states weakens and they 
become ill-defined. 

We have also modelled the LDOS images using the 
TBR model augmented with an s-wave approximation 
for the vacuum tails of the wave functions from each 
resonance state (Appendix A) . In Fig. ^3 we show the 
TBR-LDOS images of the three resonance states for the 
ad-trimer. The calculated images of the two first reso- 
nance states in Fig. [21 are rather well reproduced by this 
model. Note that the second state in the TBR model 
is anti-symmetric and its amplitudes on the atoms are 
strictly localized to the edge atoms. The absence of a 
strict nodal plane in the corresponding image is caused 
by the resonance broadening so there is still a contribu- 
tion from the first resonance state at the energy of the 
second resonance state. At the resonance energy of the 
third state, the TBR-LDOS image gives a too large con- 
tribution from the center adatom compared to the cal- 
culated LDOS image indicating that the resonance am- 
plitude on the center adatom is relatively smaller to an 
edge adatom than in the TBR model. Note that the 
LDOS images do not directly reflect the resonance am- 
plitudes on each adatom. For instance, in the TBR model 
the only difference between the resonance amplitudes on 
the adatoms between the first and the third resonance 
state is the sign of the amplitude of the center adatom 
and its magnitude is twice of that for an edge adatom. 
The large suppression of the contribution of the center 
adatom to the TBR-LDOS image at the energy of the 
third resonance compared to that of the edge adatom is 
then caused by the destructive interference amongst the 
contribution from the center adatom to the LDOS and 
those from the edge atoms in contrast to a constructive 
interference of these contribution to the LDOS at the 
energy of the first resonance state. 

For the low-energy states of longer chains with relative 
long wavelengths compared to the tip-surface distance, 
the resonance broadening results in an enhancement of 
the LDOS at the edge atoms. This effect is illustrated in 
Fig. El by the results for the TBR-LDOS profiles at an 




FIG. 10: Calculated local density of states images of a 
Au trimer on NiAl(llO) within the tight-binding, resonance 
model. First, second and third state at (a) 1.01, (b) 1.65, 
and (c) 2.29 eV, respectively. Same tip-surface distance as in 
Fig.H 



energy corresponding to the seventh (n — 7) resonance 
state of a chain with 20 atoms for three different values 
of the resonance broadening 7. At the lowest value of 7 
essentially only the n = 7 state contributes to the LDOS 
and there is only a minor enhancement of the LDOS at 
the edge atom whereas the number of contributing states 
increases with 7 and the enhancement increases at the 
edge atoms. At the largest value of 7, the nodal planes 
are no longer discernable. Alternatively, this effect may 
be understood simply as a reduction of the phase coher- 
ence length of an electron propagating along the chain 
by the resonance broadening. 



V. COMPARISON WITH EXPERIMENTS 

The unoccupied, resonances states of chains of Au 
atoms of various lengths on a NiAl(llO) surface were 
characterized by scanning tunnelling microscope both in 
an imaging and spectroscopy mode by Nilius and cowork- 
ers^. Here we are making a comparison of the calcu- 
lated resonance energies for the ad-chain Y with their 
experimental data. In addition, the "particle-in-a-box" 
(PIB) model used in the experimental analysis is clari- 
fied and justified by making a direct connection to the 
tight-binding resonance (TBR) model. 
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FIG. 11: Calculated Profiles of model local density of states 
profiles of a 20 Au atom chain on NiAl(llO) along Y within 
the tight-binding, resonance model. The profile is along the 
ad-chain with the same tip-surface distance as in Fig.0 First, 
second and third profile is at an energy of 1.39 eV correspond- 
ing to the n = 7 resonance state with a broadening 7 equal 
to (a) 0.0, (b) 0.20, and (c) 0.40 eV, respectively. The circles 
indicate the positions of the atoms in the ad-chain. 



FIG. 12: Energy variation of the lowest energy state in the 
ad-chain along the Y direction as a function of number of 
adatoms. The solid line is the result from the tight binding, 
resonance model using the effective tight binding parameters 
obtained from the calculated dispersion, whereas the dashed 
fine is the result using the "particle-in-the box" model with 
parameters obtained from the tight binding resonance model. 
The solid d iamonds are the experimental values taken from 
Refs.mH 



The calculated dispersion of the resonance states for 
the infinite, ad-chain along the Y direction is in agree- 
ment with the measured dispersion from the ad-chains 
with eleven (Aun) and twenty (AU20) adatoms. The 
dispersion for the AU20 ad-chain was obtained from an 
analysis of the characters and the energies of the observed 
resonance states in the PIB model. For the Aun ad-chain 
the second, third and fourth resonance state dominated 
the oscillations in the ^ spectra at 1.3-1.4, 1.7 and 2.5 V 
from which the number of nodal planes and correspond- 
ing wavelengths of these states and the box length could 
be determined. A fit to a free-particle dispersion gave an 
effective mass of 0.4±0.1m e . For the AU20 ad-chain, sev- 
eral resonance states contributed to the ^ spectra and 
the resonance energies, wavelengths and box length were 
extracted from a fit of the observed STS profiles to a su- 
perposition of several neighboring eigenstate densities of 
the particle in the box. The resulting effective mass of 0.5 
m e is in close agreement with the result from the Aun 
ad-chain. The calculated value of 0.65 m e compares fa- 
vorably with these measured values for the effective mass. 

Another more direct observation of the electronic res- 
onance state energies of the chain was the measurement 
of the energy, e m i n (N), of the low-energy peak in the 
STS spectra as a function of the number of atoms N 
in the ad-chain. e m in(N) was found to follow closely 
the L~ 2 dependence obtained in the PIB model where 
L is the box length. Note that the relation between L 
and N is not well-defined in the PIB model, and the 
best fit was obtained using L w (N+?)a. The calcu- 
lated value of 0.69 eV for e mi „(oo) is very close to the 
band onset of 0.68 eV obtained from an extrapolation of 
the measured values to N — > 00. For finite ad-chains, 



we can make a comparison of the calculated data with 
measured e m in(N) by using the result from the TBR 
model for e min (N) = e min (co) - 4ti sin 2 ( 2( ^ r +1) ). As 
shown in Fig. ^1 the agreement of the results from this 
model using the effective TB parameters for the infinite, 
ad-chain is in excellent agreement with the experimental 
data. Note that the difference between the results of the 
TBR and the free-particle limit is not significant. 

This good agreement between theory and experiment 
for the Au chains is somewhat fortuitous and is not 
granted by the use of unoccupied Kohn-Sham states. 
For example, in the case of the resonance states of the 
adatom, the ad-dimers and the ad-trimer the energy dif- 
ference between calculated and measured resonance en- 
ergies is typically up to 0.3 eV as shown in Fig. ^ As 
discussed in Ref. 0, a large part of this energy difference 
for the resonance in the single adatom can be accounted 
for by the Stark shift introduced by the electric field of 
the sample bias. However, in the case of the ad-trimer, 
the Stark shifts on the resonance energies are found to 
be small, that is, 0.014 to 0.03 eV for an external field of 
about 0.1 V/A. 

The success of the PIB model in the analysis of the 
observed ^ spectra and profiles, and some of its con- 
ceptual problems are clarified by the TBR model. For 
example, in the PIB model the physical meaning of the 
one-dimensional box and wave functions are not clear. 
For wave-vectors well within the LBZ boundaries, the 
TBR model gives a free-particle-like dispersion of the 
PIB model. As shown in the Appendix, the envelopes 
of the resonance state amplitudes on the adatoms in the 
TBR model are the same as the wave function in the 
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PIB model with L — (N + l)a. Furthermore, as shown 
in Fig. ^] the profile of the LDOS of a single, low- 
energy resonance state along the AU20 ad-chain in the 
TBR model is close to the density of the corresponding 
state in the PIB model. Note that for resonance states 
with higher energies and shorter wave lengths, the LDOS 
densities is enhanced at the end points. 

VI. CONCLUDING REMARKS 

We have carried out a density functional study of the 
unoccupied, resonance states in linear Au atom struc- 
tures on a NiAl(llO) surface. The primary objective has 
been to understand the character and the origin of these 
states, which were revealed and studied by scanning tun- 
nelling microscopy and spectroscopy (STS). The prime 
focus of this study has been infinite, chains of Au atoms 
along the two different orientations, [001] and [110], on 
the NiAl(llO) surface. The substantial difference in in- 
teratomic distances between these two different orienta- 
tions of the ad-chains illustrates the effects of interatomic 
interactions. From a study of the evolution of the reso- 
nance states in an ad-chain from the states of the isolated 
chain, we find that the resonance states derives from the 
6s states of the Au adatoms, which hybridize strongly 
with the substrate states and develop a substantial p-like 
polarization perpendicular to the surface. 

The calculated resonance states and their local den- 
sity of states images were analyzed in a tight-binding, 
resonance (TBR) model. This analysis provides physical 
insight about these states. The origin behind and the 
distinction between direct and substrate-mediated inter- 
actions between the resonance states are clarified in the 
TBR model. The LDOS at the end atoms of finite ad- 
chains is shown to be enhanced by two different physical 



Here cj creates an electron in an adatom state at site i 
of the adatom structure with energy £q and is the 
direct interaction (hopping) term between two adatoms 
at site i and j, V^i = V? is the interaction (hopping) 
term between an adatom state at site i and a substrate 
state /j, with energy e M and an electron in this state is 
created by c^. Note other atomic states can be included 
among the substrate states. The result for the adatom 
Green function for a single site^S generalizes directly to 
the multi-site model as, 



G tJ (e) = (^(e)- 1 ^ - tf) - E«(c)(l - Jy)) -1 (A.2) 



mechanisms. Finally, the physical meaning of "particle- 
in-box" model used in the analysis of STS for Au ad- 
chains is revealed in the TBR model. 

The calculated effective mass of free-particle-like dis- 
persion of resonance states for the infinite, Au ad-chain 
along the [001] direction is in good agreement with the 
measured effective mass, as obtained from STS of reso- 
nance states in the ad-chain with 11 and 20 Au atoms. 
The calculated energy of the band bottom is in good 
agreement with the measured energy as obtained from an 
extrapolation of the measured energies of finite ad-chains 
to an infinite ad-chain length. The principal limitations 
of describing the unoccupied resonance states in the Au 
adatom structures by unoccupied Kohn-Sham states are 
revealed by the energy differences of about 0.1- 0.3 eV 
between calculated and measured resonance energies for 
the ad-dimers and the ad-trimer. 
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APPENDIX 

The proposed tight-binding, resonance model of the 
resonance states in the linear adatom structures is based 
on a multi-site Newns- Anderson (NA) modeli2*22>. The 
Hamiltonian for this model is given by, 



(A.l) 

I 

where the Green function, C?j(e), for a single adatom at 
site i is given by, 

G i (e) = ((e-e$ ) )-X ii (e))- 1 (A.3) 

The effect of the substrate is represented by an energy- 
dependent and complex self-energy as, 

The resonance of the single Au adatom on the 
NiAl(llO) surface at the site i is now modelled as a sim- 
ple pole at cq — «7 in the Gj(e) for the 6s state with a 



J 



H = £ 4°^ + E^cjc, + h.c.) + 5]^(V,c|c, + h.c.) + £ e M ct c, 

i i,j i M M 
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pole strength Z. Note that in this case the magnitude of 
Z will be substantially less than unity. Now we assume 
that for energies around e Q + 17 the energy dependence of 
the off-diagonal part of (e) can be neglected so that 

Gij(e) w Z((e - e + i 7 )% - ty) _1 (A.5) 
where ty is a renormalized hopping term given by 

% = Z(*i?+£«(eo)) (A.6) 

This hopping term has a contribution Zty from the di- 
rect interaction between 6s atomic states renormalized 
by the resonance strength and a contribution ZHij(eo) 
from substrate-mediated interactions. Here, we restrict 
tij to nearest-neighbor interactions. Using these approx- 
imations, the multi-site adatom Green function reduces 
to an analytic form given by, 

G ij {e) = ZY J (A.7) 

where the energies ej, and the amplitudes Vfci are the 
same as in the nearest neighboring TB model for a chain 
of states and are given by, 

e k = e + 2Re(ii)cos(fca), (A. 8) 

ipki oc sin(fc|fi|) . (A. 9) 



where t\ is the nearest neighboring hopping term and 
x.[ are the positions of the N adatoms in the linear ad- 
chain with nearest neighboring adatom-adatom distance 
a, and k = (jv+i)n ' n = 1>2, ...AT. In contrast to the 
standard TB model, t\ can have a negative imaginary 
part from the substrate-mediated interactions resulting 
in a /c-dependent broadening 7^ given by, 

7 fc = 7 - 2Im(<i ) cos(fca) . (A. 10) 

Finally, note that ^st\ < will tend to increase the broad- 
ening with decreasing k and that this model contains the 
ad-dimer and the ad-trimer as special cases, correspond- 
ing to N — 2 and 3, respectively. 

Using this TBR resonance model, we have also gener- 
ated a model LDOS at the tip apex tq and an energy e 
by a superposition of localized wave functions as, 

p(f ; e) = I E M(f~ *<l t)\ 2 7r{{e _ e l ) 2 +7 2 ) 

(ATI) 

Here we have neglected any fc-dependent broadening. At 
the tip apex, we use an s-wave approximation for the 
vacuum tail of 4>(f; e) given by, 

4>(f;e) oc exp(— nr)/r, (A. 12) 

where k is the decay of the wave function in the vacuum. 
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